Perfectly conducting channel on the dark surface of weak topological insulators 
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A weak topological insulator (WTI) bears, generally, an even number of Dirac cones on its surface; 
they are susceptible of doubling, while on the surface of a certain orientation it shows no Dirac cone. 
On this "dark" surface of a WTI, we predict the existence of a single pair of isolated ID perfectly 
conducting channels that forms either a closed loop or a segment of a line. The former is associated 
typically with a single atomic-layer-thick island formed on the dark surface, while the latter is shown 
to be the consequence of a pair of crystal (screw) dislocations terminating on the dark surface. 
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The specificity of the three-dimensional (3D) topolog- 
ical insulator is that it can be weak [J-Q. Initially, this 
weak topological phase has not been paid much atten- 
tion under the pretext that it is not more than a mere 
adiabatic deformation of an ordinary insulator. Yet, it 
has later turned out that in this so-called weak phase a 
topological non-triviality is not actually weak but rather 
hidden @, 0]. In contrast to its counter part, the strong 
topological insulator (STI), which exhibits a single pro- 
tected Dirac cone, the weak topological insulator (WTI) 
bears either an even number or, depending on the way 
the crystal is cleaved, no gapless Dirac cone on its sur- 
face. And for this very reason, the WTI is often mistaken 
for possessing no robust topological character. 

The topological non-triviality ensuring the emergence 
of gapless Dirac cones on the surface of topological in- 
sulators is encoded in its non-trivial or inverted band 
structure, that is to say, stems from a non- trivial fea- 
ture in the reciprocal space. In a WTI at least one of 
the three so-called weak Z2 indices 1/1, ^3 Q remains 
non-zero, reflecting its non- trivial band structure. The 
status of the STI and WTI could be inverted when this 
reciprocal-space feature is interconnected with a topolog- 
ically non-trivial real-space character such as a disloca- 
tion of the crystalline structure in the real space. In the 
examples shown in Refs. @, a WTI exhibits a pro- 
tected gapless state bound to such a dislocation, while 
this is not necessarily the case with STIs. 

The real-space geometry, or the shape of the sample 
plays also an important role in the spectrum of the sur- 
face state 0. Indeed, except for the case of an isolated 
infinitely large single surface, the surface state is not ac- 
tually gapless in the strict sense; the gaplessness is not 
immune to finite-size gap opening. And the magnitude 
of this finite-size energy gap is qualitatively different in 
different geometries. In the simplest and most commonly 
used geometry of a slab such a finite-size energy gap is 
exponentially small, and may practically be negligible ex- 
cept in extremely thin film samples In prism-shaped 
(isomorphic to a cylinder, or a filled torus) [ToHIl or cu- 
bic (spherical) [13| STI, the magnitude of the size gap is 
much enhanced, decaying only algebraically as a function 




FIG. 1: (Color online) Upper and lower terraces separated by 
a step of height h are examples of the dark surface of WTI. 
On the lower terrace is an island formed with the same height. 



of the linear dimension of the system. 

In contrast to a STI, a WTI can have "dark" surfaces 
on which no gapless Dirac state appears. Such dark sur- 
faces can be made by cleaving the bulk crystal in the di- 
rection specified by the Miller indices (^1^2 ^3), or equiv- 
alently, by placing the surface normally to the "weak vec- 
tor" v = (^1, z/2, ^3). As a consequence of the existence of 
these dark surfaces WTI samples can exhibit size effects 
different from that of a STI of the same shape [8]. Typi- 
cally, the size effects show even/odd feature with respect 
to the number of layers stacking in the direction of v. It 
has also been pointed out that a WTI is not that weak 
either against disorder [H|, Q, [H| • Here, in this Letter we 
highlight another striking manifestation of a hidden non- 
trivial character of the WTI phase. The work is indeed 
motivated by an experimental discovery of such a WTI 
and the possibility of cleavage in the plane normal to z/, 
i.e., realizing a dark surface of WTI [161 ]. 

We start by demonstrating that a closed loop of per- 
fectly conducting channel (PCC) is formed at the pe- 
riphery of an atomic-layer-thick terrace of height corre- 
sponding to an odd number of atomic layers formed on 
a dark surface of WTI. Experimental observation of such 
atomic-scale terraces formed on an otherwise flat surface 
of a topological insulator by STM measurements have 
been also reported previously [Tt( ] . 

Let us consider the following standard Wilson-Dirac 
type effective Hamiltonian for 3D topological insulators, 
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FIG. 2: (Color online) A pseudo-gapless ID mode emergent 
along the periphery of an island. The island has a size of 
£x X £ y = 9 X 11, and an altitude h = 1. The systems's sizes 
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16, N z = 7. 



here implemented on the cubic lattice as 

#buik = ATxd^smk^ +r z m(k), 
m(k) — mo + 27712/1(1 — cos fc M ), 



(i) 



where two types of Pauli matrices a and r represent 
real and orbital spins, and /i = x,y,z. The parameter 
A represents the strength of normal hopping, and deter- 
mines the aperture of the surface Dirac cone. To realize a 
WTI with a specific choice of weak indices v = (0, 0, 1), 
we choose the mass parameters to be anisotropic. In 
the following demonstrations they are chosen such that 
mo = —1, rri2 X = m2y = 0.5 and rri2 Z = 0.1 in units 
of the hopping parameter A [8]. In the actual simula- 
tions on-site random potential is also introduced as in 
Kef. [15j ]. The global shape of the sample is also as- 
sumed to be tetragonal apart from terraces and islands; 
see below, i.e., the underlying lattice is, unless otherwise 
mentioned, terminated in all the three (x-, y- and z-) 
directions, and in principle could accommodate a sur- 
face state on all of its tetragonal faces, but because of 
the weak vector pointing to the ^-direction that excludes 
Dirac cones from the surfaces normal to this direction, 
the system exhibits pseudo-gapless surface states only on 
the four side surfaces. 

Let us now consider a situation in which the top sur- 
face is not completely flat, showing either atomic-scale 
terrace-like step structures or an isolated island-shaped 
region as depicted in Fig. [H When the surface is normal 
to z/, the vertical wall of this step structure is parallel 
to z/, on which one can hypothesize formation of gapless 
2D surface states (two Dirac cones). Then, of course, 
contributions from these two Dirac cones must be super- 
posed to cope with the boundary condition that the wave 
function is not allowed to leak onto the dark terraces [i| . 

Fig. [2] and Fig. [3] (a) demonstrate that there appears 
indeed a pair of ID modes almost localized at the periph- 
ery of such a terrace- or an island-shaped region, when its 
height h (measured in units of the lattice constant) is odd 
(h = 1, 3, 5, • • • ). These ID modes are, in contrast to 2D 
surface states, immune to backward scattering as a con- 
sequence of their helical nature, and realize an example 
of what is often called "perfectly conducting channels" 
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FIG. 3: (Color online) Profile of the lowest energy wave func- 
tion in atomic-scale terraces of a WTI. (a) When h = 1, ID 
modes appear at the edge of the terrace, (b) when h = 2, 
the wave function has a vanishing amplitude along the edge 
of the terraces. N x = N y = 16, N z = 7. 



[18|. The wave functions associated with these ID chan- 
nels shown in Fig. [2] and Fig. [3] (a) are barely influenced 
by the presence of weak disorder; here, the strength of 
disorder is chosen as W/A = 1 in the notation of [sj]. 
The above numerical simulations have been performed 
in systems of a rather simple shape of the island and of 
the terrace. While, the appearance of a ID PCC is con- 
sidered to be generic, and equally applicable to systems 
that are smoothly connected to the present one as far as 
h is odd. 

If h is even, the situation is much different as a result of 
the even/odd feature of finite-size gap generation charac- 
teristic to the surface of WTI @. This even/odd feature 
is associated with the fact that a WTI is sometimes de- 
scribed as stacked layers of 2D quantum spin Hall states, 
each exhibiting a pair of helical edge modes. These he- 
lical edge modes, when stacked to form a helical surface 
state, tend to couple into gapped pairs, leaving only a 
single gapless channel when the number of stacked layers 
is odd. A more quantitative k • p-type argument devel- 
oped in Ref. [8[ suggests when applied to the terrace- or 
island-shaped region that the pseudo-zero-energy states 
acquire a finite energy, E± = ^ 2 (/^+i) 7r ^ wnen its height 
h is even. Here, the pseudo-gapless ID helical modes are 
completely kicked out of the low-energy spectrum [see 
Fig. [3] (b)]. The case of a two-story (h = 2) terrace 
is an example in which this (parity of) /^-dependent size 
effect is most accentuated, giving rise to a huge finite- 
size correction to the energy of pseudo-zero modes as, 
E± = ±A 7r, i.e., on the order of 1 in units of A. 

The protected ID modes emergent when h is odd, typi- 



FIG. 4: (Color online) A numerical experiment: verification 
of the perfectly conducting character of the ID channels cir- 
culating the periphery of the island (case of h = 1). The 
spectrum of the ID modes is plotted against the inserted flux 



FIG. 5: (Color online) A finite segment of ID gapless modes 
penetrating into the bulk WTI at the ends of the segment 
along a screw dislocation line with b — bz (b — dbl). The 
two dislocation lines are separated by a distance of ( = 8. 
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cally one, are not entirely immune to finite-size gap open- 
ing, but of a different, less destructive type. On one hand, 
these counter propagating modes are always in the op- 
posite spin state (spin-to- momentum locking), and this 
very fact makes them immune to backward scattering by 
the impurities. While, on the other hand, the quanti- 
zation axis of these spins can and does indeed rotate as 
an electron propagates along this circular ID orbital, in 
such a way that this quantization axis is normal to the 
tangent of this ID orbital (spin-to-surface locking); a sit- 
uation analogous to the one occurring at a cylindrical 
surface of STI To|- 1_2] . As a result, while an electron re- 
volves once around a circular orbit, its spin also performs 
a complete 2tt rotation, altering the orbital part of the 
boundary condition from periodic to anti-periodic. This 
leads, tracing exactly the same path as in the case of the 
half-odd integral quantization of the orbital angular mo- 
mentum L 7 in a cylindrical STI 0, EH, to a finite-size 
correction E± = ±^7r, where A is a constant of order A, 
and s = 4£ is the circumference of the island. 

To double check that the pseudo-gapless ID modes in 
the case of h being odd realize indeed a single pair of per- 
fectly conducting channels, and that they are immune to 
backward scattering by the impurity, we have performed 
the following numerical experiment. In a situation as de- 
picted in Fig. 2, in which a ID ring-shaped distribution 
of the wave function is realized along the periphery of an 
island, we gradually introduce a magnetic flux tube (of 
magnitude ^>) threading the island such that ideally the 
flux would not touch the electron. The boundary con- 
dition on the side surfaces are here chosen to be doubly 
periodic, i.e., a slab geometry is employed to avoid 2D 
states on the side surfaces. Fig. [4] depicts the result of 
such a numerical experiment in which the spectrum of 
the ID states is plotted against <j). Here, </> = 27r^- is 
an Aharonov-Bohm phase associated with the external 
flux measured in units of the flux quantum <l>o- One 
can clearly observe that (i) the spectrum is doubly de- 
generate at (j) = and at <j) = 7r, (ii) as (j) is introduced, 



say at <\> = 0, this two- fold degeneracy is lifted in such 
a way that one state is shifted upward, and the other 
downward, indicating the existence of counter propagat- 
ing channels. Note also that (iii) at <fi = and at <p = 7r, 
the pairs are swapped, and (iv) the (^-dependence of each 
spectrum is almost perfectly linear. These features are 
perfectly in accordance with our hypothesis that a pair 
of ID helical modes are formed at the periphery of the 
island. 

The spectrum of the ID helical modes in the presence 
of a flux can be written most generally as 



s 



2~£ +AE > 



(2) 



where AE represents a global shift of the spectrum as- 
sociated with the breaking of the particle-hole symme- 
try due to on-site random potentials. At (j) = 0, the 
two-fold degeneracy of the spectrum occurs such that 
E^~\o) = £?ij +1 (0), while at <fi = 7r, after these partners 
are shifted in the opposite directions, perfectly linearly 
to </>, each state gets paired again but with a different 
partner such that E^~\tt) = e[tJ(tt), i.e., in such a way 
that the relative ordering of the partners is shifted by one 

We have seen so far examples in which pseudo-gapless 
protected ID modes appear on a dark surface of WTI, 
forming either a closed loop, e.g., at the periphery of an 
island, or a line, e.g., at the boundary of two terraces, 
crossing the dark surface completely, only to merge with 
a 2D side-surface state at both ends. In the remainder of 
this Letter we focus on another type of ID modes emer- 
gent on the dark surface, typically as a finite segment 
of a line of length (. We demonstrate that this type of 
ID modes are associated with the presence of dislocation 
lines in the crystal, and indeed the two ends of the seg- 
ment coincide with the termination of a dislocation line 
at the dark surface (see Fig. [5]). 

To convince oneself that this is certainly the case let 
us first turn around the argument in the following way. 
Recall that a pair of gapless ID helical modes are in- 
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FIG. 6: (Color online) The "flux experiment" applied to a 
configuration with a pair of screw dislocation lines. 

duced along a dislocation line introduced in the bulk WTI 
0, > when the Burgers vector b specifying the disloca- 
tion satisfies the relation, b • M = tt mod 27r, where 
the vector M is defined in terms of the weak indices 
and the three reciprocal lattice vectors Gi, G2, G3 as 
M = \{y\G\ + V2G2 + V3G3). The question to be ad- 
dressed here can then be reformulated in the following 
way. What happens to the single pair of gapless ID he- 
lical modes that emerges along a dislocation line, when 
the dislocation is terminated at a dark surface of WTI? 
In any realistic setup it is natural to assume that a dis- 
location line is terminated at the surface of the sample 
at both ends, unless it forms a closed loop. Here, we 
consider a dislocation line terminated by the two dark 
surfaces, such as the ones on the top and on the bottom 
(see Fig. [5]). Since these dark surfaces are normal to 
z/, the condition for the existence of gapless ID helical 
modes is satisfied, only when the introduced dislocation 
is a screw dislocation with a dislocation line parallel to b. 
When such a screw dislocation is terminated at the sur- 
face of a sample, it necessarily accompanies a mismatch 
of atomic layers by a magnitude of b = |6|, forming typi- 
cally a step structure, analogous to the one that appears 
between the two terraces. 

Let us come back to the original point of view in which 
our focus is on a finite segment of the ID helical modes 
emergent on the dark surface. Notice that a single pair 
of such Dirac states cannot be confined, in principle, in 
a finite segment (Klein tunneling). The only way (out) 
for such a segment of ID modes to appear on the dark 
surface is that it pair-annihilates with another ID helical 
modes along a dislocation line that terminates at the dark 
surface. This is explicitly demonstrated in Fig. [5l 

The last question to be addressed is whether these 
continuous ID channels partly hosted by a dislocation 
line, combined with a segment of ID modes on the dark 
surface realize also a perfectly conducting channel. For 
this purpose we have performed a numerical experiment, 
analogous to the one for the island geometry, but here, 
the flux tube is inserted between the two dislocation lines. 
In Fig. [6] one can clearly recognize the four characteristic 
features supporting the formation of a pair of ID helical 



modes observed in the case of the island. Here, in the 
present geometry this confirms that different parts of the 
helical modes, one on the dark surface, and the other 
along a dislocation line, are indeed smoothly connected 
to form a ID perfectly conducting channel. Here, s in 
Eq. ([2]) is given by s = 2(£ + N z ), the circumference of 
the ID channel. 

A single pair of pseudo-gapless ID helical modes also 
appears along a dislocation line inserted in a STL How- 
ever, when the pair is hosted by a STI, it necessarily 
terminates in the "sea" of the single-cone Dirac state at 
the surface where surface sensitive measurements are ap- 
plicable. Here, we have provided a more experimentally 
feasible situation to detect a dislocation solely from in- 
formation on the surface, by predicting that termination 
of a ID perfectly conducting channel on the dark surface 
of WTI accompanies a screw dislocation line in the bulk 
host crystal that terminates at the dark surface. 
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